NAS A/TM— 2006-2 14359 



Three-Dimensional Field Solutions for Multi-Pole 
Cylindrical Halbach Arrays in an Axial Orientation 


William K. Thompson 

Glenn Research Center, Cleveland, Ohio 


September 2006 



NASA STI Program ... in Profile 


Since its founding, NASA has been dedicated to the 
advancement of aeronautics and space science. The 
NASA Scientific and Technical Information (STI) 
program plays a key part in helping NASA maintain 
this important role. 

The NASA STI Program operates under the auspices 
of the Agency Chief Information Officer. It collects, 
organizes, provides for archiving, and disseminates 
NASA’s STI. The NASA STI program provides access 
to the NASA Aeronautics and Space Database and its 
public interface, the NASA Technical Reports Server, 
thus providing one of the largest collections of 
aeronautical and space science STI in the world. 
Results are published in both non-NASA channels and 
by NASA in the NASA STI Report Series, which 
includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant phase 
of research that present the results of NASA 
programs and include extensive data or theoretical 
analysis. Includes compilations of significant 
scientific and technical data and information 
deemed to be of continuing reference value. 
NASA counterpart of peer-reviewed formal 
professional papers but has less stringent 
limitations on manuscript length and extent of 
graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies that 
contain minimal annotation. Does not contain 
extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, often 
concerned with subjects having substantial 
public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific and 
technical material pertinent to NASA’s mission. 

Specialized services also include creating custom 

thesauri, building customized databases, organizing 

and publishing research results. 

For more information about the NASA STI 

program, see the following: 

• Access the NASA STI program home page at 
http://www.sti.nasa.gov 

• E-mail your question via the Internet to 
help@sti. nasa.gov 

• Fax your question to the NASA STI Help Desk 
at 301-621-0134 

• Telephone the NASA STI Help Desk at 
301-621-0390 

• Write to: 

NASA STI Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076-1320 



NASA/TM— 2006-214359 



Three-Dimensional Field Solutions for Multi-Pole 
Cylindrical Halbach Arrays in an Axial Orientation 


William K. Thompson 

Glenn Research Center, Cleveland, Ohio 


National Aeronautics and 
Space Administration 


Glenn Research Center 
Cleveland, Ohio 44135 


September 2006 



Acknowledgments 


The author gratefully acknowledges Dennis J. Eichenberg, Jeffrey Juergens, and Dawn C. Emerson, NASA Glenn Research 
Center, for assistance with the development of the analytical models; Christopher A. Gallo, NASA Glenn, for contributions to 
the graphics; and Mark Christinin, Ansoft, Inc., for assistance with the finite-element model. 


Trade names and trademarks are used in this report for identification 
only. Their usage does not constitute an official endorsement, 
either expressed or implied, by the National Aeronautics and 
Space Administration. 


This work was sponsored by the Fundamental Aeronautics Program 
at the NASA Glenn Research Center. 


Level of Review. This material has been technically reviewed by technical management. 


Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076-1320 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22161 


Available electronically at http://gltrs.grc.nasa.gov 



Three-Dimensional Field Solutions for Multi-Pole Cylindrical 
Halbach Arrays in an Axial Orientation 

William K. Thompson 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

1. Introduction 

Uses for the so-called Halbach array of permanent magnets have grown in number in recent years. 

The salient feature of the Halbach array is the unique placement of individual permanent magnets such 
that the B field is concentrated on one side of the array and canceled on the other. This useful and 
intuitively efficient property exists for both linear and cylindrical arrays (ref. 1). In addition to the 
inventor’s original designs for particle beam focusing mechanisms and undulators, one may now find 
Halbach arrays in a number of applications, including high-performance motors and generators (ref. 2), 
frictionless passive magnetic bearings and couplers (refs. 2 and 6) and magnetically levitated trains (ref. 3). 

This article presents three-dimensional B field solutions for the cylindrical Halbach array in an axial 
orientation. This arrangement has applications in the design of axial motors and passive axial magnetic 
bearings and couplers. The analytical model described here assumes ideal magnets with fixed and 
uniform magnetization. The model also assumes a sufficiently large number of magnets (N m > 16) 
comprise the Halbach array so that the angular span of each individual magnet is kept small. This permits 
modeling its magnetization as arising from a sum of four surface currents. The field component functions 
are expressed as sums of 2-D definite integrals that are easily computed by a number of mathematical 
analysis software packages. The solutions are found to be sinusoidal functions of angular position (with 
additional harmonics present at axial distances that are small compared to the magnet thickness), 
exponential functions of axial distance from the magnets and more complex functions of radial position 
that must be computed numerically. The analysis is verified with sample calculations and the results are 
compared to equivalent results from traditional finite-element analysis (FEA). The field solutions are then 
approximated for use in flux linkage and induced EMF calculations in nearby stator windings by 
expressing the field variance with angular displacement as pure sinusoidal function whose amplitude 
depends on radial and axial position. The primary advantage of numerical implementation of the 
analytical approach presented in the article is that it lends itself more readily to parametric analysis and 
design tradeoffs than traditional FEA models. 

2. Magnetic Field Theory of The Axial Halbach Array 

Figure 1 shows the cylindrical Halbach array in an axial orientation. The term “axial Halbach array” 
will be used from this point forward to refer to this configuration. The array depicted is comprised of 
N m = 32 sector shaped permanent magnets with inner radius /•/, outer radius r 2 , and axial thickness T. 
Assume that N m > 1 6 so that the angle in radians spanned by each magnet is small compared with 2 n. 
There are four magnets per Halbach wavelength in the angular {(j)) direction. Each sector in the array has 
an index s, where s = [0, 1 , ... N m - 1 ]. The magnets each have a magnetization M = ±B,//j 0 whose direction 
is indicated by the arrows. B r is the remanent magnetization of the permanent magnet material. 

Field calculations require the definition of two overlapping coordinate systems, one Cartesian and the 
other cylindrical. The cylindrical r</> plane aligns with the Cartesian xy plane, so that the z coordinates of 
each system are identical. The +x axis of the Cartesian system lies at ^ = 0 in the cylindrical system. The 
field solutions will be found at any arbitrary point in space, which defines a vector in the cylindrical 
system given by r = (r r + (fn/) + z z ). The vector giving the location of the integration variable is 
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Figure 1. — The axial Halbach array for N m = 32 magnets. Arrows indicate the direction of magnetization 
for each individual magnet. This particular arrangement will concentrate the B field below the ring and 
cancel it above the ring. 


designated by primes, i.e., r '= (r'r + <f>' <j> + z'z ). The distance between these two points is given by 
Green’s function 


2 + r ,2 + (z -z ') 2 - 2 rr ' cos(^ - ^ ') ^ 

The array lies in the Cartesian system such that the +z axis corresponds to the axis of rotation and the 
flat bottom faces of the sectors lie in the xy plane. Hence, zj = 0 and r? = T. The 5 = 0 magnet is selected 
to be magnetized in the axial (+z) direction and the +x axis bisects the flat bottom surface of this magnet. 
Magnets lying on the -x, +y, and —y axes also bear this same direction of magnetization. The 5 th magnet 
has a central radial axis which lies at angle </> = /3 S = 2 re s/N,„ Each magnet subtends an angle in the 
(f) direction of /T - fit = 2 tt I N m radians, where /T and /?/ are the locations of the side faces. 



The magnetic field component solutions are expressed in cylindrical coordinates as 

B{r,4,z) = {B r ,B t ,B,} 

and these may be transformed to Cartesian coordinates using 

B(x,y,z) = {B r cos^ - sm<p, B r sm<j) + B^ cos(p, B _ } 


( 2 ) 


( 3 ) 
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Using the principle of linear superposition, one may express the aggregate field components as the 
sum of individual contributions from each sector in the array. For any axially magnetized (±z) magnet in 
the array, the B field components may be calculated from the vector potential and expressed as sums of 
definite integrals in two of the spatial dimensions as determined previously (ref. 4): 
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Note that +M is used for +z magnetization and -M is used for -z magnetization. In a similar manner, 
expressions for the field components for an individual transversely magnetized sector may be expressed 
in the following equations, which are derived in the appendix of this article. 
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Here the assumption of sufficiently large N m becomes particularly important as the magnetization is 
not truly azimuthal but linear and perpendicular to the radial axis of the magnet at (p - J3 S . Note that +M is 
used for +(f) magnetization and -M is used for -</> magnetization. 

The contributions of each magnet in the array add to give the aggregate solution. Starting at the 5 = 0 
magnet and traveling around the array in the +</> direction, the magnetization directions repeat the pattern 
{ + z ,+ 0 ,-z,-0 , + z ...}. 

Therefore, the magnetizations for each sector in the array are 

Axial Case: M = {-l)"Mz (6a) 

Transverse Case: M = (-1 )" M(j> (6b) 

A new indexing variable n has been defined to account for the interleaving of axially and transversely 
magnetized sectors. One may express the angular position of each axially magnetized sector as 


/?» = 


and the left and right faces of the sector correspond to angles at 


o (4 n - l)/r o {An + 1 )n 

P\ — 5 P 2 — 

N ' N 


Similarly, each transversely magnetized sector is located at an angular position 


PM)- 


(An + 2)7: 


and the left and right faces of the sector correspond to angles 


o {An + 1 )n o (4« + 7>)k 

P — ? P 2 — 

N N 


this yields the field components for the entire collection as 
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These expressions define the B field components at any point in space. They may be easily 
implemented in a variety of commercial mathematical analysis software packages. We used Mathematica 
v5.2 (Wolfram Research, Inc). Special care must be taken to observe the signs of each of the terms, which 
are dictated by the j variable of summation. 


NASA/TM— 2006-214359 


5 



3. Results and Validation of Field Solutions 


The analytical expressions for the B field components given in section 2 have been coded in 
Mathematica. This product permits easy numerical implementation of the definite integrals using the 
built-in NIntegrate[] function. The working precision for the numerical integrations was set to 50 digits 
and the accuracy goal to five decimal places. Equivalent FEA models were also developed in Maxwell 3D 
vlO (Ansoft, Inc.). The validation method compares the results of these two independent models of the 
same axial Flalbach array. 

The selected design parameters for the simulations are: 

r, = 1.0” (25.4 mm), r 2 = 2.0” (50.8 mm) 

T = z 2 - zi = 0.25” (6.4 mm) 

N m = 32 magnets, 4 magnets per Halbach wavelength 
B r = HoM= 1.5 T (NdFeB-55 rare earth permanent magnets) 

These parameters match those of an axial magnetic bearing model currently under development at 
NASA Glenn Research Center (ref. 5). 

Figures 2(a), (b), and (c) compare the radial, angular and axial dependence, respectively, of the B, 
component of the analytical and FEA models at an axial gap distance z = -0.05 in. (—1.3 mm). 

Radial dependence of B r at z = -0.05 in. 


0.6 

0.4 

0.2 
K 

c 0.0 
CQ 

- 0.2 
-0.4 
- 0.6 

0 012345678 

cj), integer multiples of rinlN m 



Angular dependence of B r at z = -0.05 in. 



Axial dependence of B r at <j) = 4n/N m 



z, in. 

Figure 2. — Axial Halbach array at selected radial and axial distances and angular locations. Results of numerical 
integration of the analytical expressions are shown as smoothed lines. Equivalent FEA results are shown as 
markers at selected points, (a) B r versus r. (b) B r versus <|>. (c) B r versus z. 
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Figures 3(a), (b), and (c) compare the radial, angular and axial dependence of the B 0 component of 
the analytical and FEA models at an axial gap distance of z = -0.05 in. (—1.3 mm). Figures 4(a), (b), and 
(c) compare the radial, angular and axial dependence of the B : component of the analytical and FEA 
models at an axial gap distance ofz = -0.05 in. (—1.3 mm). Figures 5(a), (b), and (c) plot the three B field 
components in the first quadrant of the xy plane at an axial gap distance ofz = -0.05 in. (—1.3 mm). 
Finally, Figure 6 compares the azimuthal dependence of the three B field components at a larger gap 
distance of z = -0. 15 in. (-3.8 mm) and at the radial location which produces the maximum field: 
r = 1.35 in. (34 mm). 


Radial dependence of at z = -0.05 in. 



Angular dependence of at z = -0.05 in. 



<j), integer multiples of nnlN m 


Axial dependence of at <]) = 1 . 5jr/A/ m 



z, in. 

Figure 3. — Axial Halbach array at selected radial and axial distances and angular locations. Results of numerical 
integration of the analytical expressions are shown as smoothed lines. Equivalent FEA results are shown as 
markers at selected points, (a) B r versus r. (b) B r versus c|>. (c) B r versus z. 
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Figure 4. — Axial Halbach array at selected radial and axial distances and angular locations. Results of numerical 
integration of the analytical expressions are shown as smoothed lines. Equivalent FEA results are shown as 
markers at selected points, (a) B z versus r. (b) B z versus <|>. (c) B z versus z. 



NASA/TM— 2006-214359 








Figure 5. — 3D plots of FEA results for B field components B r ( a), Ba (b) and B z (c) of the axial Halbach array in the 
first quadrant of the plane at an axial distance of z = -0.03 in. below the magnets. 




(b) Radians 


Figure 6. — 2D plots of Br, Bq and Bz versus angular position (<|>) at r = 1.35 in. (34 mm) and z = -0.15 in. (4 mm) as 
computed by FEA (a) and analytical (b) models. Both models indicate that at sufficiently large values of gap (g) the 
field components exhibit nearly perfect sinusoidal behavior with negligible harmonic components. 
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4. Discussion 


The plots indicate agreement to within 1 0 percent of peak values between the numerical 
implementation of the analytical model presented in this article and an equivalent FEA model. The radial 
component of the field B, is the weakest of the three and has significant magnitude only at the two radial 
edges of the array. For most practical applications, this field component is of no consequence. 

Both the B,/, and B z components have higher amplitude and greater spatial extent and achieve their 
maximum amplitudes at radial distances r = 1.35 in. (34 mm) and r = 1.317 in. (33.5 mm), respectively. 
The B : component is flatter than B,/, over the radial extent of the magnets, as can be seen by comparing 
Figures 4(a) and 5(a). Flowever, both B,f, and B z must be treated as functions of radial position, B/r) and 
B-(r) when computing flux linkage and induced emf and current in nearby stator windings. 

At gap values which are small compared to the magnet thickness, T, all three field components 
contain significant harmonics versus angular position (j). The angular frequency of the fundamental is a>i = 
8 n/N m as can be seen from figures 3(b), 4(b), and 5(b). At an axial position ofz = -0.05 in. (-1.3 mm) 
distortion of the sinusoid arises mostly from a 5th harmonic. As z — > 0 one finds the presence of even 
higher order harmonics. Flowever, for most practical applications the field is used to compute flux linkage 
in a winding of relatively large spatial extent where the gap is sufficiently large such that all harmonic 
content may be neglected. For larger values of z, only the fundamental sinusoid of frequency f remains as 
shown in figure 7. Flere the gap value is -0.15 in. (-3.8 mm). 

Post has determined that the axial dependence of the field components is an exponential function 
(ref. 3) of the form B = B a e k ~ where k = 2n/X and A is the Flalbach wavelength, i.e., the width of four 
magnets along the direction of travel. For the axial Flalbach array, the direction of travel is azimuthal and 
k is not a constant as in (ref. 3) but a function of radial position given by 

k(r)=NJ4r (11) 

Summarizing all of this, we can write field component equations in the manner described by Post 
(ref. 3) using the form 


= B (/)0 sin (f) e~ kz 

(12a) 

= B :o cos (j) e~ k 

(12b) 


with the understanding that the B,i, t B :o and k are all functions of r. B ^ and B :o give the field component 
strengths of the fundamental at the surface of the Flalbach magnet array (z = 0) and at radial location r. 
This gives 


B r/) ( r , (j), z ) = B r/io (r) sin (j) e k(r): (13a) 

B z {r,(j>,z) = B zo (r) cos </>e~ k(r)z (13b) 

as the working equations for the axial Flalbach array. From these, dynamic analysis of a rotating array 
may follow by substituting (j) = where (/>, is the initial angular position of the array and co is the 

mechanical frequency of rotation. The analysis is simplified if the array is initially positioned with $• = 0. 
Dynamic analysis is explored further in a related article (ref. 5). 

Furthermore, simplifying assumptions may be made regarding the effect of relatively distant magnets 
from the spatial location of interest. The field effects of magnets more than two Flalbach wavelengths 
away from a spatial location near the array may be considered negligible. Factoring this into the 
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programming of the numerical computation of the analytical model may reduce computation times 
drastically, especially for very large N m . Finally, in a related article (ref. 6) Post simplifies further by 
assuming an average value over the span of a nearby stator winding for the peak field, rather than the 
function of radial position indicated above. This produces a closed form approximation to the field that 
may be sufficiently accurate for a variety of applications. 

5. Conclusion 

This article presents analytical expressions for the B field solutions for an axial Flalbach array of 
permanent magnets. The analytical expressions are easily implemented in numerical analysis software 
packages. Validation of the analytical model by finite element analysis shows agreement between the two 
methods within 1 0 percent of the peak value of the field. 
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Appendix 

Derivation oiB Field Solutions for Transversely Magnetized Sector Magnets 


A sector of permanent magnetic material with inner radius r t and outer radius r 2 is positioned such 
that the center axis of the bottom surface coincides with an azimuthal angle (</> = J3 S ) and the bottom face 
of the sector lies in the xy plane (z/= 0). Although only field effects from a single sector are considered 
here, the sector is known to be part of a complete circular array of N m identically-sized magnets and 
therefore subtends an angle P=2ji/N m . The angular component <j> is defined over the range [0, 2 n\. 

We proceed with the derivation of B via the vector potential A in a similar method as Furlani 
(ref. 4). Figure 7 shows the geometry of the magnetization in terms of the spatial integration variables 
r (j)' and z . The magnetization vector M lies in the transverse plane. For a sufficiently large number of 
magnets in the total array, the angle subtended by each individual magnet is small. In this case the 

magnetization may be approximated as having a (j) component only, given by 


M - ±M(j> 


(14) 


since J = V x M and the magnetization has no net circulation, there is no volume current density. The 
magnetization therefore arises from surface current densities on the top, bottom, inner and outer surfaces 
of the magnet sector, where 


j M = M x « 


(15) 



Figure 7. — Coordinate systems for the derivation of the component B field 
expressions for a single transversely magnetized sector at an angular 
position, |3 S . M points in the direction of $ at the centerline of the sector, 
where (|) = p s . Primed quantities indicate variables of integration. 
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The four surfaces are described analytically as 


top surface = < 


bottom surface = < 


inner surface = < 


outer surface = < 


r < r 1 < r 
1 2 

A <</>'< fi 2 

Z' = z 2 

r < r ? < r 
1 2 

z' = z t 
r' = r 

^<r<j3 2 

Zj < Z ' < Z 2 

r ' = r 2 

^<<t>'<j3 2 

z, < z ' < z. 


The four surface normals are given by 


n = < 


+z (top) 

—z ( bottom ) 
-f (inner) 
+r (outer) 


and the four surface current densities are given by 


Mr 

(top) 

—Mr 

(bottom) 

Mz 

(inner) 

—Mz 

(outer) 


We will derive B from the vector potential function, A, using 

B =Vx2 


and Green’s function 


G(r,r') ■ 


1 


V 


r" + r '"+ (z — z ')“ - 2rr 'cos(^ - 


(16a) 


(16b) 


(16c) 


(16d) 


(17) 


(18) 


(19) 


( 20 ) 


which gives the distance between the spatial location of interest and the variables of integration. The 
vector potential arises from volume and surface current densities, but the volume density has already been 
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shown to be zero. Therefore the vector potential will arise from four surface integrals of surface current 
densities given by 


A(r) = ^~ [ — da' 
An ' I r - r '] 


( 21 ) 


Substituting equation (16) into this expression for j M gives 


4<o=±— £<-iy 


4 n 


7=1 


r i Pi 

H 

A 


-r'dfdr' 


Zl Pi 

u 

, z i A 


-r'dfdz' 


( 22 ) 


but r is a function of the angular position variable <f>'. Substituting the transformation 


r — cos (j)' x + sin^’y 


(23) 


and arranging terms gives 


A s (r) = ± 


An 



r 2 Pi , , - 

f 

+ 

Z '=Z ; 

r 2 Pi . , , - 

J 

J J \r-r'\ 

~ 

L'-. pi 1 1 

z'= Zj 

z 2 Pi 


f f = r'dtb' dz' 

J J 1 r - r ' 1 


L z i Pi' 1 

- r'=Vj 


(24) 


which is the vector potential at the point of interest in terms of the Cartesian unit vectors. Projecting these 
terms back into cylindrical coordinates using 


x = cos <j)r - sin (fxp (25a) 

y = sin </>r + cos (jxj) (25b) 

one may collect like terms and simplify using angle sum and difference trigonometric identities to obtain 
the cylindrical coordinate components of the vector potential function as 


A rJr) = ± 


An 


!<-» 


>2 pi 


7+1 


— cos(0 ~ <j>') , 


7=1 


JJ 7T vr ,,r ' r'dfdr' 


>\ Pi 


r - r 


J Z'=Z: 


(26a) 
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(26b) 


4, ( F)=±^i(-ir 


4tt j=\ 
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z 2 P 2 
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J J I r - r ' I 
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z i A 


(26c) 


and the final step is to take the curl of A in cylindrical coordinates to obtain the individual components of 
the B field. 


15 8 _ P 0 M ^, 
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